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Abstract. We present the formulation of the problem of the coherent dynamics 
of quantum mechanical two-level systems in the adiabatic region in terms of the 
differential geometry of plane curves. We show that there is a natural plane curve 
corresponding to the Hamiltonian of the system for which the geometrical quantities 
have a simple physical interpretation. In particular, the curvature of the curve has the 
role of the nonadiabatic coupling. 


1. Introduction 

Two-level quantum systems (TLS), being the simplest of nonsimple quantum systems 
[1], can be used to study a wealth of phenomena in quantum physics. Basically, it 
describes systems with two discrete states, for example the spin degrees of freedom of 
spin-1 particle, but is often used to describe systems only effectively, as in quantum 
optics where it serves as a model for resonant excitation of an atom [2]. More recently, 
with the advent of the held of quantum information, understanding the dynamics of 
TLS that act as the quantum information carriers and are therefore dubbed qubits, has 
become even more important in order to both understand the fundamental physics and 
gain control over the information processing tasks [3]. 

Considering modern physics more generally, the appreciation of the fact that 
the physical phenomena should not depend on the particular coordinate description 
employed by a physicist has become one of the central principles and has lead to 
increasingly refined and abstract geometrical tools and language with which to express 
the physical laws [1]. In contrast, in this paper we shall discuss the dynamics of a 
driven TLS in a more elementary differential geometrical setting, but nevertheless try 
to exploit fully the geometrical character of the resulting theory. 

Coherent driving of the system with classical external fields and parameters leads 
to time-dependencies in the matrix elements of the system Hamiltonian. We can take 
the Hamiltonian to be real symmetric and traceless without any loss in generality [5]. 
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This means that a particular model, i.e. external driving in question, is described in 
general by two different functions of time. These are usually called as the detuning 
and the Rabi frequency, stemming from the quantum optical setting, but they can have 
any physical origin, such as the two independent spatial components of the external 
magnetic field in the case of a spin-1 particle. The models where the evolution of 
the state vector is exactly solvable are rare, despite the long history of the problem 
[a m m e], From a mathematical point of view, a TLS Hamiltonian, or the detuning 
and Rabi frequency defining it, can be thought of as defining a plane curve with time 
variable as its parameter. We show that there is a natural way of defining the plane 
curve from these functions in a way that the relations between the resulting geometrical 
quantities and the basic physical variables conventionally associated with the TLSs 
exhibit certain simplicity. Despite the elementary character of the results, these are not 
usually explicitly discussed in treatments of quantum dynamics of TLS. There are many 
interesting papers mmm in a similar vein that relate the Hamiltonian of the TLS 
and the geometry of plane or space curves but the way the relation is done differs from 
our approach. 

The paper is organised as follows. In the section [ 2 ] we first give the mathematical 
formalism and definitions used to describe coherently driven TLS in 12. II and connect the 
theory with the differential geometry of plane curves in 12.21 We give some examples of 
the usefulness of this connection in subsection 12.31 The discussion in [3] will conclude this 
presentation. The basic mathematical results of plane curves are given as an Appendix. 

2. Mathematical Formalism 

2.1. Two-level Quantum Systems 

Generally, the coherent time evolution of a quantum system is given by the Schrodinger 
equation, in units where h — 1, 

i d t if(t) = H(t)i/>(t), (1) 

where H ( t ) is the Hamiltonian operator with time-dependent components, acting on the 
Hilbert space of the system. For TLS, the space consists of two-dimensional complex 
state vectors ijjft) = (c+(i), c_(t)) T , where c + {t) and c_(t) are the complex probability 
amplitudes of the states formed by the natural basis of the vector space. Physically 
these basis states, denoted by tp±, usually refer to the two system state vectors in the 
absence of interaction and the basis is usually called bare or diabatic. In this basis the 
Hamiltonian is denoted by H d and takes the form 

H d (t) = H(t ) • a (2) 

(m m \ 

- [ m —a (t)) ’ 

where the first line contains the shorthand form involving the field vector H(t) = 
(f2(t), 0, A(i)) and the Pauli matrices. The function A(t) is called detuning and f 1(t) 
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is called the Rabi frequency. This is the most general form needed for TLS, as any 
Hermitian two-by-two matrix can be transformed into this form by suitably redefining 
the phases of the basis vectors [5j. Instead of the held vector given in the cartesian 
coordinates, it is also useful to consider the Hamiltonian in the spherical coordinates, 
in which case 


H d (t) = pit) 


cos ( 0 (f)) sin( 0 (£)) 

sin( 0 (f)) —cos (0(t)) 


(3) 


where p{t) = \JA 2 (f) + fi 2 (t) is the length of the held vector and 9 it) is the angle it 
makes with the z-axis, so that Q/A = tan(0). The eigenvalues of the Hamiltonian, the 
quasi-energies, are given by 


E±(t) = ±p(t)- 


(4) 


The corresponding eigenstates of the Hamiltonian (j3J) and (l3|) , which can be chosen 
to be real, read 


X+(t) = ± 



X-{t) = ± 



(5) 


The adiabatic basis, which is formed by these time-dependent eigenvectors (choosing 
sign to be the same for both eigenvectors), is obtained by x± — Up>± where the time- 
dependent unitary transformation is given by 

cos(-rp') -sin(^) \ 

sin(^p) cos(^jl) J 

The Schrodinger equation in the adiabatic basis reads 

. d_ f a+(t ) \ = / p[t) i lit) \ ( a+it) 
d t y a_(f) J y -i 7 (t) -pit) J y a-(t) 

where if = a + y + + a_x_ and the fact that the transformation is time-dependent induces 
a gauge term in the adiabatic Hamiltonian that couples the adiabatic basis states. This 
term is called adiabatic coupling 7 (f) and is given by 




lit) = “(x+|X-(t)) 

_ A(t)(l(t) - A(t)Q(t) 

2 (A 2 (t) + H 2 (f)) 


where the overhead dot stands for time derivation. When the Hamiltonian changes 
slowly, it is well known that the transitions between the adiabatic states vanish and the 
evolution is called adiabatic. We take a_ (tj) = 1 as the initial condition and define the 
final transition probability as P = \a + (t/) | 2 where fj and tf refer to initial and final 
times, respectively. One can often take = — 00 and tf = +00 but we consider also 
finite initial and final times. The condition for adiabatic evolution is given by [2] 

|A(«)f 2 (() - A(«)fi(«)l « [A 2 (() + tf(()] 3/2 . 


(9) 
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We consider here the generic case where the functions A(i) and Q(t) have different 
zeros, so there is only avoided crossings in adiabatic energies. However, the diabatic 
energy levels, given by ±A(£), can cross, a fact that is exploited in the rapid adiabatic 
passage techniques [13]. A physically relevant assumption is that the Rabi frequency 
that couples the diabatic states vanishes when t —> ±oo. Then, or more generally, when 
the condition 


-C |A(i)|, t —* tij, 


( 10 ) 


is satisfied, the adiabatic and diabatic basis coincide at the initial and final times. 
However, one should note that depending on the character of the detuning function, the 
diabatic and adiabatic basis states may swap labels in relation to each other. That is, 
it can happen that X-iff) — P+ although x~(U) — and in this case the transition 
from the state to state (p + can be obtained with an adiabatic time evolution. In 
spherical coordinate description, this is obtained when 

9(t f ) - d(ti) = 2 J 7(*)d t = 2tt (n + 0 , (11) 


where n is an integer and the bases coincide in the initial time when we choose 9(tf) = 0. 
Other interesting cases then are 9(tf) = 2im and 9(tf) — n (n + |) which correspond 
to complete population return and formation of equal superposition, respectively, in the 
adiabatic evolution. 

It should be also noted that the assumption of the non-crossing adiabatic levels 
allows us to define the change of the time variable, 

s{t) = f p(x)dx, (12) 


Jo 

which means that the time is measured in terms of the dynamical phase that is 
accumulated. This transformation of the time variable proves to be very important 
in what follows, and it has the effect of normalizing the held vector in ([3]) and the 
eigenenergies (J4j), which become E±(s) = ±1. Therefore, we could choose the detuning 
and the Rabi frequency as 


A(s) = cos(6 ) (s)) (13) 

H(s) = sin(6 ) (s)), (14) 


which shows that there is only one function needed to describe any model of coherent 
TLS dynamics, namely 9(s). 


2.2. Connection to Plane Curves 

A curve considered as a connected set of points in plane is a very intuitive concept 
and one can easily picture it in ones mind. Therefore, it is worthwhile to connect 
the driven TLS problem of the previous subsection with some suitably defined plane 
curves. To this end, it is convenient to consider explicitly a parametrized plane curve 
a{t) = (x(t),y(t)). The basic results from the theory of differential geometry of plane 
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curves that are needed in the following are given, along with the further references, in 
the Appendix. We choose to associate with each diabatic Hamiltonian (J3D a plane curve 
whose components are given by 

x(t) = / A (u)du, y(t ) = f Q(u)du, (15) 

Jo Jo 

the parameter t being the physical time and usually we take t £ (—oo,+oo). Note 
that the so-called Hamiltonian curves defined by Berry in [12J and used to study driven 
quantum systems are the velocity, or the (unnormalized) tangent, vectors to the curves 
defined by (fT5l) . With our definition, the curvature of the curve is simply related to the 
functions in the Hamiltonian and given by 

27(f) 


n[a](t) = 


p(f) ’ 


(16) 


as is obvious when we look at the formulas (J8]) and (1A.4D . The adiabatic coupling 
is directly proportional to curvature and the adiabatic condition © can be simply 
translated to 

k[ck](s) 1. (17) 

Furthermore, the speed of the curve is now just given by p{t) and the existence 
of only avoided crossings means that the curve (1151) is regular and its unit-speed 
parametrization is given by the change of the time variable (fT2l) . In this parametrization, 
the curvature and its relation to the adiabatic coupling simplifies further as the 
denominators in (fT6l) and (1A.4I) become unity. Then from (1T4l) it follows that curvature 
is given by 

d 9(s) 


n[a](t) = 


ds 


(18) 


which shows that the angle in (TT4|) is the same as the so-called turning angle of the 
curve (g5D , up to a constant (see figured]). 


2.3. Applications and examples 

The curves that are formed from the known soluble models are typically very simple and 
are not bounded in any finite region of the plane. This is due to fact that we usually have 
U — —tf — —00 while the adiabatic coupling, and therefore the corresponding curvature, 
differ appreciably from zero only near avoided crossings (which are usually chosen to 
happen near t — 0) and so the curves have well-defined lines as their asymptotes. Figure 
|T] depicts the situation for Landau-Zener [1] [7] E] E| and parabolic [Il][I5yT6j models. It 
also shows whether or not the model swaps the labels between the diabatic and adiabatic 
states, as discussed in connection with (TTTT) . 

However, we can take the opposite point of view and consider some well-known plane 
curves with interesting geometric properties from the beginning and then study what 
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kind of Hamiltonians they give rise to. For example, the ellipse is given by parametric 
equations 

An 

xi(t) = —sin (ut) (19) 

UJ 

yi(t) = — (1 - cos (cut)), (20) 

UJ 

where A 0 , Ho and uj are positive constants. The way the ellipse is defined here places 
it completely on the positive half of the y-plane and the respective magnitudes of the 
constants Ao and Ho determine the minor and major axes, i.e., whether the ellipse is 
squeezed in X\- or yi-direction. The corresponding time-dependent Hamiltonian can be 
described by the held vector in the (e 1; e 2 , e 3 )-coordinate system 

Hi(t ) = (fi 0 sin(cuf), 0, A 0 cos(ujt)) T , (21) 


which also gives an ellipse rotating in ei-e s plane. Although the angular frequency uj 
can be a physical parameter, it does not affect the geometrical properties, so we set it 
equal to unity. In particular, one should note that the value of uj does not affect the 
adiabaticity (see equation (1A.4I1 ). We furthermore restrict our time to the one-period 
interval t G [0, 27r]. This diabatic Hamiltonian is then an example of zero-pulse model 
with two level crossings at t — | and t = -y- with any parameter values of A 0 and fi 0 . 
The adiabatic and diabatic bases coincide at initial and final times. This can, of course, 
be obtained directly from (l2Tj) but a nice way is also to imagine the ellipse ([20]) and to 
note the connection between the function 0(s), defining the Hamiltonian in (TT4l) . and 
the turning angle, which makes it immediately obvious that 0(0) = 0(27r) = 0. The 
curvature of the original curve (1201) is given by 

A 0 fi 0 


Ki(t) = 


(Aq cos 2 (f) + fig sin 2 (f)) 


3 / 2 ’ 


( 22 ) 


a) Landau-Zener model b) Parabolic model 




x(t)/A 0 x(t)/A 0 

Figure 1 . Examples of plane curves defined by TLS model Hamiltonians when 
t G (—3,3): a) Landau-Zener model, A (t) = A 0 t, fl(t) = fi 0 , b) The double-crossing 
parabolic model, A(i) = Aq (t 2 — l), f l(t) = fig where Ag and fig are constants. 
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Pi{t) = \J Aq cos 2 (f) + sin 2 (f). (23) 

The latter, of course, gives the eigenenergies of the Hamiltonian and also the adiabatic 
coupling is directly read from these expressions via The corresponding quantities 
are plotted in figure [2l 


Energy levels 



Coupling 




3 4 5 6 


0 1 


2 3 4 

t 





Figure 2. Examples of the diabatic (full line) and the adiabatic (dashed) energy level 
structure of the Hamiltonian defined by the ellipse (12(1 are plotted on the left panel 
while the corresponding couplings are depicted on the right panel. In the uppermost 
case we have Ao > Ho and Ho = 0.2, in middle panel we have the case of a circle with 
radius Ho = 0.5 and in the lowermost picture we have Ho > Ao, Ho = 1.5. In all of 
the plots A 0 = 0.5. 


The basic connection made between the plane curves and Hamiltonians in 12.21 that 
was based on a simple observation about the similarities between the mathematical 
expressions of different quantities, may seem rather trivial. The point, however, is that 
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now one can make use of the many strong differential geometric results for plane curves. 
One of the most celebrated results, of global character, is the four-vertex theorem (FVT) 
Id which says that the curvature of a simple closed curve (a circle exluded) has to have 
in total at least four points of local maxima and minima. Because of the association of 
the nonadiabatic coupling as curvature and the fact that multiple peaks in the coupling 
can lead to interference effects because of the distinct times when the transition happens, 
it is interesting to study the implications of FVT, namely, that for certain type of models 
the distinct peaks in the adiabatic coupling can not be made to go away. 

As a simple plane curve, the ellipse obeys the FVT and it has exactly four vertices, 
at t — 0, |, 7T, ^y. Which ones correspond to maximas and which minimas, depends on 
the parameters A 0 and V2 0 , but in any case k(0) = k(tt) and ac(|) = «(^y) and these 
values are 



(24) 


(25) 


The adiabatic limit of this model is not immediately obvious. If O 0 > A 0 , the maxima of 
curvature are at points t = 0, it and minima at t — |, ^y. These maximum peaks become 
higher and more narrow when the ratio fio/Ao gets bigger. On the other hand, keeping 
Ho fixed, these maxima of curvature become smaller when we increase the value of Ao- 
However, this has the effect of increasing the value of the curvature at the minimum 
points. Indeed, when A 0 = Ho? the curvature is a constant function with non-zero value 
A^ 1 as the curve corresponds to a circle of radius Ao- The adiabatic limit for the TLS 
defined by this circle is obtained as its radius tends to infinity. Also for a general ellipse 
we see that the adiabatic condition (fI71) is reached, for example, if we increase both fi 0 
and Ao but keep the ratio fio/Ao constant. The case when we further increase the value 
of A q , so that A 0 > O 0 , is basically similar to the case O 0 > A 0 , but now the maxima 
and minima of the curvature have swapped places. This is particularly interesting when 
we compare the evolution in diabatic and adiabatic bases. 

In diabatic basis, the basic structure of the energy levels and corresponding 
couplings remains the same in all parameter regions. It has two level-crossings at fixed 
times when the coupling is also maximal, while the coupling is zero in the initial and 
final times, see figure EJ In adiabatic basis, we have to consider different parameter 
regions. Although the bases coincide at initial and final times in all parameter regions, 
the structure of the adiabatic levels and couplings depend on the respective magnitudes 
of A 0 and fi 0 . When A 0 > Ho the avoided crossings of the levels happen at the same 
instant as the crossings in the diabatic basis and also the adiabatic coupling is maximal 
there. So in both bases the transitions are located near these two points. When Ho > Ao, 
the maxima of the curvature swap with the minima as discussed above, and there is only 
one avoided crossing and a peak in the adiabatic coupling between t t and tf at t — n. 
However, the adiabatic coupling obtains its maximum also at the initial and final times. 

As another example we can consider the limagon curve which is not simple so that 
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FVT is not valid in this case. It is given in the parametric form by 

X 2 (t) = ( acos(ut ) + b ) cos(aA) (26) 

t/ 2 (t') = (acos(aA) +b)sin(ujt), (27) 

where t G [0, 27t] and we again have three parameters a, b and oj, all appearing now in 
both components. The latter parameter is set to unity for the same reason as in the 
previous example. Of course there are special curves defined by (|27|) which correspond 
to special parameter values. For example, when either a = 0 or b = 0 we have a circle 
and the considerations of the adiabatic limit in the previous example apply. On the 
other hand, if a = b, the equations ()27l) define a cardioid, which has a singularity in the 
curvature function at t — n, and so we cannot expect to obtain an adiabatic evolution 
over the whole curve. The curve (|27l) gives rise to a Hamiltonian with the detuning 
function and Rabi frequency as 

A 2 (t) = —bsin(t) — asin(2f) (28) 

0 2 (t) = bcos(t ) + acos(2f). (29) 

Curvature and speed in the case are given by 

2 a 2 + b 2 + 3afecos(f) 


K 2 (t) = 


(30) 


(a 2 + b 2 + 2 ab cos(f)) 3 / 2 ’ 
p 2 (f) = \J a 2 + b 2 + 2a6 cos ( t ). (31) 

We consider the case where a > b > 0, so that the curve is not simple. It is still regular 
and has an inner loop, see figure |3] The FVT does not apply now and curvature has 
only one maximum, at t — n and 

1 


k 2 (tt) = 


1 + 


(32) 


a — b \ a — b , 
and the adiabatic limit is obtained when \a — b\ tends to infinity. 

We have chosen for convenience the initial time as t* = 0 but of course any one- 
period time interval, i.e. T = 2n, would suffice to trace the limagon curve. With the 
current convention, the diabatic levels start at degeneracy and with maximal diabatic 
coupling. A similar discrepancy between the description in different bases that was met 
in the first example is also present here (see figure |3|) . In adiabatic basis, there is only 
one avoided crossing, at t — tt, and the corresponding peak in the coupling between 
the adiabatic basis states. This coupling is smallest in the initial and final times when 
the basis states are equally weighted superpositions of the diabatic basis states because 
d(0) = |. Meanwhile, the diabatic levels go through three crossings during the evolution. 

As a final example, we note that a Hamiltonian of the form 

A 3 (f) = acos(nt + 5) (33) 

^3 (f) = 6cos(f), (34) 

where a, b , n and 6 are real parameters gives rise to the Lissajous curve 


x 3 (t) = — sin {nt + S ) 


n 


y 3 (t) = f)siii(f). 


(35) 

(36) 
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Again, this curve has well-known properties and it can be in general quite complicated. 
The simple special cases include the circle (n — 1, a = b, 5 = §), a line (n — 1, 5 = 0) 
and a parabola (n = 2, 5 = |). In any case it is bounded around the origin by a box 
with sides — and 2b and when n is rational, n — j, it is closed and the corresponding 
Hamiltonian is periodic. The curve is nevertheless highly sensitive to the changes in the 
values of n and the nominator k and denominator l give the number of ” lobes” in vertical 
and horizontal directions, respectively. So arbitrarily small deviation of n from unity 
will affect greatly to the Lissajous figure compared to the simple elliptical shape when 
n = l, see figure [H This means that small deviations in the angle frequencies of the 
different trigonometric functions in the Hamiltonian (1HT]1 give very different behaviors 
for the quantum system, given long enough time, as one can associate each ” lobe” with 
some maximum of the curvature. 




0.0 0.5 1.0 1.5 2.0 2.5 0 1 2 3 4 5 6 



0 1 2 3 4 5 6 


x 


t 


Figure 3. Left: Limagon curve of equation (1271) with a = 1.5, 6=1. Middle: Energy 
levels corresponding to Hamiltonian (l2f)l) . Diabatic levels are depicted in full line while 
the adiabatic levels are dashed. Right: The diabatic (full line) and the adiabatic 
(dashed) coupling of (l29l) . 


Lissajous curve 




! .5 f V V V V V V V V V 
0 10 20 30 40 50 60 


Coupling 


\pJt 

A A A A AW 

w 

t jimMfvv v 

' 


6 1 -- — .— 

0 10 20 30 40 50 60 


x 


t 


Figure 4. Left: An example of Lissajous curve of equation (l36l) . Middle: Energy 
levels corresponding to Hamiltonian (l34l) . Diabatic levels are depicted with a full blue 
line while the adiabatic levels are purple and dashed. Right: The diabatic (full line, 
blue) and the adiabatic (dashed, purple) coupling of (1291) . The parameters in all of 
the plots are a = b = 1.5, S = ^ and n = 0.9. It should be noted that had the n been 
equal to unity, these plots would be similar to the middle panel of figure [2] 
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We have discussed the general formulation of the dynamics of the time-dependent two- 
level quantum systems and pointed out how there is a natural correspondence in the 
theory of plane curves. It is remarkable that by simply defining components of a curve 
parametrically as integrals of the basic functions of Hamiltonian, namely detuning of 
the energy levels and coupling between the basis states, we get simple geometrical 
interpretations for the basic physical variables. Taking this point of view, adiabatic 
coupling is curvature, eigenenergy is the speed of the curve and the polar angle of the 
Hamiltonian in (EJ is the turning angle of the curve, for example. Also the coherent 
adiabatic dynamics is completely given by the geometry of the curve in a clear and 
intuitive manner. Non-adiabatic transitions happen near the points where the curve 
bends most and crossings of the eigenenergies will induce singularities in the adiabatic 
coupling. To consider dynamics outside adiabatic region one can, for example, consider 
curves defined analogously in the higher-order superadiabatic bases of p3] , where the nth 
and (n+ l)th order bases take the role of the diabatic and adiabatic bases, respectively. 

One important consequence of this formulation is that now one can apply the 
mathematical results of the theory of plane curves to the dynamics of two-level systems 
in a direct way. The fundamental theorem of plane curves gives a constructive way 
to obtain interesting Hamiltonians with given adiabatic couplings. Also, we discussed 
how the four-vertex theorem gives conditions for the approach of the adiabatic limit in 
the case of Hamiltonians which correspond to simple closed curves. Unlike the curve of 
the parabolic model in figure |T] which in the adiabatic limit tends to a line, the total 
curvature of a closed plane curve is always a multiple of 2ti so it can tend to the adiabatic 
limit somewhat differently, by smoothing the peaks in the curvature but also increasing 
in length. 

Many of the parametrizations of the closed plane curves involve trigonometric 
functions and we considered three different periodic level-crossing models motivated by 
the differential geometry of such plane curves. Each of these exhibit interesting features, 
for example, when it comes to parameter dependence and to the level structures and 
couplings when compared in different bases. These models were not directly obtained 
from a particular physical problem but this is not a drawback since nowadays one can 
modify the time-dependencies of the external fields at will, with the advent of new laser 
technologies, for example. It should be also noted that many simple models can be 
obtained for certain special parameter values for the studied models. 
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We introduce here just the basic definitions and properties of plane curves that suffice 
to our purposes of relating them with time-dependent TLS as discussed in the main 
text. The definitions and notations here are mostly standard and more information can 
be found, for example, from [f?, H2J ED]- 

A parametrized curve is defined as a (piecewise-) differentiable function a : I —> R n , 
where I = (a, b ) is an open interval in R, either finite or infinite. Furthermore, we often 
differentiate the curves freely so in effect we actually assume the curves to be smooth. 
As we are interested in curves in the plane, we set n = 2. The plane curve can be 
considered as a parametrized vector given in a Cartesian coordinates as 


a{t) = (x(t),y(t)), 


(A.l) 


where x(t) and yft ) are two real functions and t G I. The derivative of a curve is obtained 
naturally by componentwise differentiation and a and a are called the velocity and 
acceleration of a curve a, respectively. The norm of the velocity vector, v(t) = ||b(f)||, 
is called speed. For a curve to be well-behaved it is required that v(t) 0 for all t and 
such curves are called regular. 

Two distinct functions can trace the same point set on the plane, so it may not 
always be immediately clear when two curves are actually the same, i.e., they differ only 
by the parametrization. If we have two curves, a : (a, b) — > R 2 and [3 : (c, d) —> R 2 and 
there exists a differentiable function h : (c, d ) —¥ (a, b ) such that h(t) > 0 (h(t) < 0) for 
all c < t < d and (3 = a o h, we say that f3 is a positive (negative) reparametrization 
of a. The different signs of the reparametrization are related only to the direction the 
curve is traversed. From an intuitive geometrical point of view it is also clear that any 
purely geometrical quantity associated with curves should not depend nontrivially on 
the parametrization. 

One such basic geometric quantity is the length L[a\ of a curve a : ( a,b ) —>■ R" , 
defined by 

L[a] = f \\a'(u)\\du. (A.2) 

J a 

A closely related quantity is the arc-length function of curve. Fix a number c G (a, b ) 
and let the upper limit of integration be the variable t, 


l[a,c](t) — J ||d' , ('u)||d'u. (A.3) 

For any regular curve a{t) there exists a reparametrized unit-speed curve /3(s), 
meaning that ||/3 / (s)|| = 1. The unit-speed parameter, which is essentially unique 
(up to change of origin and sign), is hereafter denoted by s. Because of the property 
/[/3, c](s) = s — c, unit-speed curves are said to be parametrized by the arc length. The 
unit-speed parametrization is very useful as many of the formulas simplify when it is 
used. 
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The most important quantity one can associate with a plane curve is its curvature 
k [a] (t). It basically measures the way the plane curve differs from a straight line, being 
identically zero only for a line and constant if and only if the curve is an arc of a circle. 
A formula for a regular curve is given by 


k [a] (t) = 


x(t)y{t) - x(t)y(t ) 


(A.4) 


(i 2 <«) + f(t)f 2 ' 

To see more clearly the meaning of the curvature function, one can associate two 
orthonormal vectors on each point of a unit-speed curve, the tangent vector t(s) and 
a vector obtained by rotating this by 7t/2, namely the normal vector n(s). Then the 
curve can be obtained from the Frenet equation 


A ( *(*) ^ = ( 0 «(*) ( t(8) \ 

ds n(s) J y — k(s) Oyl n(s) J ’ 


The content of the fundamental theorem on plane curves is that curvature 
determines the plane curve essentially uniquely, meaning up to Euclidean motions and 
reparametrizations. A curve with any desired curvature can be realized with the unit- 
speed construction 


a(s) = ( / cos 9(s)ds + c, / sin 9(s)ds + d 


9(s ) = / n(s)ds + 9q, 


(A.6) 


where c, d and 9q are constants. The function 9(s) is also called the turning angle of the 
curve and is the angle between the x-axis and the tangent vector of the curve at point 
t. 

Finally, the four-vertex theorem states that the curvature function of any simple 
closed plane curve, other than a circle, must have at least four vertices, that is, points 
where the curvature is locally minimal or maximal. One should also note that if we 
relax the condition of simplicity, the curve can have less than four vertices. 
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